Insights into the genetics of menopausal vasomotor symptoms: genome-wide analyses of routinely-collected primary care health records.

BACKGROUND: Vasomotor symptoms (VMS) can often significantly impact women's quality of life at menopause. In vivo studies have shown that increased neurokinin B (NKB) / neurokinin 3 receptor (NK3R) signalling contributes to VMS, with previous genetic studies implicating the TACR3 gene locus that encodes NK3R. Large-scale genomic analyses offer the possibility of biological insights but few such studies have collected data on VMS, while proxy phenotypes such as hormone replacement therapy (HRT) use are likely to be affected by changes in clinical practice. We investigated the genetic basis of VMS by analysing routinely-collected health records. METHODS: We performed a GWAS of VMS derived from linked primary-care records and cross-sectional self-reported HRT use in up to 153,152 women from UK Biobank, a population-based cohort. In a subset of this cohort (n = 39,356), we analysed exome-sequencing data to test the association with VMS of rare deleterious genetic variants. Finally, we used Mendelian randomisation analysis to investigate the reasons for HRT use over time. RESULTS: Our GWAS of health-records derived VMS identified a genetic signal near TACR3 associated with a lower risk of VMS (OR=0.76 (95% CI 0.72,0.80) per A allele, P=3.7x10-27), which was consistent with previous studies, validating this approach. Conditional analyses demonstrated independence of genetic signals for puberty timing and VMS at the TACR3 locus, including a rare variant predicted to reduce functional NK3R levels that was associated with later menarche (P = 5 × 10-9) but showed no association with VMS (P = 0.6). Younger menopause age was causally-associated with greater HRT use before 2002 but not after. CONCLUSIONS: We provide support for TACR3 in the genetic basis of VMS but unexpectedly find that rare genomic variants predicted to lower NK3R levels did not modify VMS, despite the proven efficacy of NK3R antagonists. Using genomics we demonstrate changes in genetic associations with HRT use over time, arising from a change in clinical practice since the early 2000s, which is likely to reflect a switch from preventing post-menopausal complications in women with earlier menopause to primarily treating VMS. Our study demonstrates that integrating routinely-collected primary care health records and genomic data offers great potential for exploring the genetic basis of symptoms.

Menopause age, reproductive span and hormone therapy duration predict the volume of medial temporal lobe brain structures in postmenopausal women.

Medial temporal lobe (MTL) atrophy is correlated with risk and severity of Alzheimer disease (AD) pathology and cognitive decline. Increasing evidence suggest that oestrogens affect the aging of MTL structures. Here we investigate the relationship between reproductive hormone exposure, polygenic scores for AD risk and oestradiol concentration, MTL anatomy and cognitive performance in postmenopausal women. To this end, we used data from 10,924 female participants in the UK Biobank from whom brain MRI and genetic data were available. We fitted linear regression models to test whether the volume of structures comprising the MTL were predicted by a) timing related to menopause, b) the use and timing of hormone replacement therapy (HRT) and c) polygenic scores for AD risk and oestradiol concentration. Results showed that longer use of HRT was associated with larger parahippocampal volumes (2.53 mm3/year, p = 0.042). A later age of natural menopause, and a longer reproductive span, was associated with larger hippocampal (6.08 and 5.72 mm3/year, p = 0.0006 and 0.0005), parahippocampal (4.17 mm3 and 4.19 mm3/year, p = 0.00006 and 0.00001), amygdala (2.10 and 2.22 mm3/year, p = 0.028 and 0.01) and perirhinal cortical (2.56 and 2.95 mm3/year, p = 0.028 and 0.008) volumes. Superior prospective memory performance was associated with later age at natural menopause, and a longer reproductive span (ß = 0.05 and 0.05 respectively, p = 0.019 and 0.019). Polygenic scores for AD risk and for oestradiol concentration were not associated with MTL volume and did not interact with menopause-related factors to affect MTL structure. Our results suggest that HRT use did not have any detrimental effects on cognition or brain structure, whilst greater exposure to reproductive hormones across time is associated both with slightly larger volumes of specific MTL structures and marginally superior memory performance, independent of genetic risk for AD and genetic predisposition for higher oestradiol levels. However, the clinical utility of maintenance of oestrogens post-menopause for brain health and protection against cognitive decline is curtailed by the small effect sizes observed.

Penetrance of pathogenic genetic variants associated with premature ovarian insufficiency.

Premature ovarian insufficiency (POI) affects 1% of women and is a leading cause of infertility. It is often considered to be a monogenic disorder, with pathogenic variants in ~100 genes described in the literature. We sought to systematically evaluate the penetrance of variants in these genes using exome sequence data in 104,733 women from the UK Biobank, 2,231 (1.14%) of whom reported at natural menopause under the age of 40 years. We found limited evidence to support any previously reported autosomal dominant effect. For nearly all heterozygous effects on previously reported POI genes, we ruled out even modest penetrance, with 99.9% (13,699 out of 13,708) of all protein-truncating variants found in reproductively healthy women. We found evidence of haploinsufficiency effects in several genes, including TWNK (1.54 years earlier menopause, P = 1.59 × 10-6) and SOHLH2 (3.48 years earlier menopause, P = 1.03 × 10-4). Collectively, our results suggest that, for the vast majority of women, POI is not caused by autosomal dominant variants either in genes previously reported or currently evaluated in clinical diagnostic panels. Our findings, plus previous studies, suggest that most POI cases are likely oligogenic or polygenic in nature, which has important implications for future clinical genetic studies, and genetic counseling for families affected by POI.

Genetically proxied therapeutic prolyl-hydroxylase inhibition and cardiovascular risk.

Prolyl hydroxylase (PHD) inhibitors are in clinical development for anaemia in chronic kidney disease. Epidemiological studies have reported conflicting results regarding safety of long-term therapeutic haemoglobin (Hgb) rises through PHD inhibition on risk of cardiovascular disease. Genetic variation in genes encoding PHDs can be used as partial proxies to investigate the potential effects of long-term Hgb rises. We used Mendelian randomization to investigate the effect of long-term Hgb level rises through genetically proxied PHD inhibition on coronary artery disease (CAD: 60 801 cases; 123 504 controls), myocardial infarction (MI: 42 561 cases; 123 504 controls) or stroke (40 585 cases; 406 111 controls). To further characterize long-term effects of Hgb level rises, we performed a phenome-wide association study (PheWAS) in up to 451 099 UK Biobank individuals. Genetically proxied therapeutic PHD inhibition, equivalent to a 1.00 g/dl increase in Hgb levels, was not associated (at P < 0.05) with increased odds of CAD; odd ratio (OR) [95% confidence intervals (CI)] = 1.06 (0.84, 1.35), MI [OR (95% CI) = 1.02 (0.79, 1.33)] or stroke [OR (95% CI) = 0.91 (0.66, 1.24)]. PheWAS revealed associations with blood related phenotypes consistent with EGLN's role, relevant kidney- and liver-related biomarkers like estimated glomerular filtration rate and microalbuminuria, and non-alcoholic fatty liver disease (Bonferroni-adjusted P < 5.42E-05) but these were not clinically meaningful. These findings suggest that long-term alterations in Hgb through PHD inhibition are unlikely to substantially increase cardiovascular disease risk; using large disease genome-wide association study data, we could exclude ORs of 1.35 for cardiovascular risk with a 1.00 g/dl increase in Hgb.

Applying a genetic risk score for prostate cancer to men with lower urinary tract symptoms in primary care to predict prostate cancer diagnosis: a cohort study in the UK Biobank.

BACKGROUND: Prostate cancer is highly heritable, with >250 common variants associated in genome-wide association studies. It commonly presents with non-specific lower urinary tract symptoms that are frequently associated with benign conditions. METHODS: Cohort study using UK Biobank data linked to primary care records. Participants were men with a record showing a general practice consultation for a lower urinary tract symptom. The outcome measure was prostate cancer diagnosis within 2 years of consultation. The predictor was a genetic risk score of 269 genetic variants for prostate cancer. RESULTS: A genetic risk score (GRS) is associated with prostate cancer in symptomatic men (OR per SD increase = 2.12 [1.86-2.41] P = 3.5e-30). An integrated risk model including age and GRS applied to symptomatic men predicted prostate cancer (AUC 0.768 [0.739-0.796]). Prostate cancer incidence was 8.1% (6.7-9.7) in the highest risk quintile. In the lowest quintile, prostate cancer incidence was <1%. CONCLUSIONS: This study is the first to apply GRS in primary care to improve the triage of symptomatic patients. Men with the lowest genetic risk of developing prostate cancer could safely avoid invasive investigation, whilst those identified with the greatest risk could be fast-tracked for further investigation. These results show that a GRS has potential application to improve the diagnostic pathway of symptomatic patients in primary care.

Detection and characterization of male sex chromosome abnormalities in the UK Biobank study.

PURPOSE: The study aimed to systematically ascertain male sex chromosome abnormalities, 47,XXY (Klinefelter syndrome [KS]) and 47,XYY, and characterize their risks of adverse health outcomes. METHODS: We analyzed genotyping array or exome sequence data in 207,067 men of European ancestry aged 40 to 70 years from the UK Biobank and related these to extensive routine health record data. RESULTS: Only 49 of 213 (23%) of men whom we identified with KS and only 1 of 143 (0.7%) with 47,XYY had a diagnosis of abnormal karyotype on their medical records or self-report. We observed expected associations for KS with reproductive dysfunction (late puberty: risk ratio [RR] = 2.7; childlessness: RR = 4.2; testosterone concentration: RR = -3.8 nmol/L, all P < 2 × 10-8), whereas XYY men appeared to have normal reproductive function. Despite this difference, we identified several higher disease risks shared across both KS and 47,XYY, including type 2 diabetes (RR = 3.0 and 2.6, respectively), venous thrombosis (RR = 6.4 and 7.4, respectively), pulmonary embolism (RR = 3.3 and 3.7, respectively), and chronic obstructive pulmonary disease (RR = 4.4 and 4.6, respectively) (all P < 7 × 10-6). CONCLUSION: KS and 47,XYY were mostly unrecognized but conferred substantially higher risks for metabolic, vascular, and respiratory diseases, which were only partially explained by higher levels of body mass index, deprivation, and smoking.

Identifying molecular mediators of the relationship between body mass index and endometrial cancer risk: a Mendelian randomization analysis.

BACKGROUND: Endometrial cancer is the most common gynaecological cancer in high-income countries. Elevated body mass index (BMI) is an established modifiable risk factor for this condition and is estimated to confer a larger effect on endometrial cancer risk than any other cancer site. However, the molecular mechanisms underpinning this association remain unclear. We used Mendelian randomization (MR) to evaluate the causal role of 14 molecular risk factors (hormonal, metabolic and inflammatory markers) in endometrial cancer risk. We then evaluated and quantified the potential mediating role of these molecular traits in the relationship between BMI and endometrial cancer using multivariable MR. METHODS: Genetic instruments to proxy 14 molecular risk factors and BMI were constructed by identifying single-nucleotide polymorphisms (SNPs) reliably associated (P < 5.0 × 10-8) with each respective risk factor in previous genome-wide association studies (GWAS). Summary statistics for the association of these SNPs with overall and subtype-specific endometrial cancer risk (12,906 cases and 108,979 controls) were obtained from a GWAS meta-analysis of the Endometrial Cancer Association Consortium (ECAC), Epidemiology of Endometrial Cancer Consortium (E2C2) and UK Biobank. SNPs were combined into multi-allelic models and odds ratios (ORs) and 95% confidence intervals (95% CIs) were generated using inverse-variance weighted random-effects models. The mediating roles of the molecular risk factors in the relationship between BMI and endometrial cancer were then estimated using multivariable MR. RESULTS: In MR analyses, there was strong evidence that BMI (OR per standard deviation (SD) increase 1.88, 95% CI 1.69 to 2.09, P = 3.87 × 10-31), total testosterone (OR per inverse-normal transformed nmol/L increase 1.64, 95% CI 1.43 to 1.88, P = 1.71 × 10-12), bioavailable testosterone (OR per natural log transformed nmol/L increase: 1.46, 95% CI 1.29 to 1.65, P = 3.48 × 10-9), fasting insulin (OR per natural log transformed pmol/L increase: 3.93, 95% CI 2.29 to 6.74, P = 7.18 × 10-7) and sex hormone-binding globulin (SHBG, OR per inverse-normal transformed nmol/L increase 0.71, 95% CI 0.59 to 0.85, P = 2.07 × 10-4) had a causal effect on endometrial cancer risk. Additionally, there was suggestive evidence that total serum cholesterol (OR per mg/dL increase 0.90, 95% CI 0.81 to 1.00, P = 4.01 × 10-2) had an effect on endometrial cancer risk. In mediation analysis, we found evidence for a mediating role of fasting insulin (19% total effect mediated, 95% CI 5 to 34%, P = 9.17 × 10-3), bioavailable testosterone (15% mediated, 95% CI 10 to 20%, P = 1.43 × 10-8) and SHBG (7% mediated, 95% CI 1 to 12%, P = 1.81 × 10-2) in the relationship between BMI and endometrial cancer risk. CONCLUSIONS: Our comprehensive MR analysis provides insight into potential causal mechanisms linking BMI with endometrial cancer risk and suggests targeting of insulinemic and hormonal traits as a potential strategy for the prevention of endometrial cancer.

Genome-wide association study meta-analysis identifies three novel loci for circulating anti-Müllerian hormone levels in women.

STUDY QUESTION: Can additional genetic variants for circulating anti-Müllerian hormone (AMH) levels be identified through a genome-wide association study (GWAS) meta-analysis including a large sample of premenopausal women? SUMMARY ANSWER: We identified four loci associated with AMH levels at P < 5 × 10-8: the previously reported MCM8 locus and three novel signals in or near AMH, TEX41 and CDCA7. WHAT IS KNOWN ALREADY: AMH is expressed by antral stage ovarian follicles in women, and variation in age-specific circulating AMH levels has been associated with disease outcomes. However, the physiological mechanisms underlying these AMH-disease associations are largely unknown. STUDY DESIGN, SIZE, DURATION: We performed a GWAS meta-analysis in which we combined summary statistics of a previous AMH GWAS with GWAS data from 3705 additional women from three different cohorts. PARTICIPANTS/MATERIALS, SETTING, METHODS: In total, we included data from 7049 premenopausal female participants of European ancestry. The median age of study participants ranged from 15.3 to 48 years across cohorts. Circulating AMH levels were measured in either serum or plasma samples using different ELISA assays. Study-specific analyses were adjusted for age at blood collection and population stratification, and summary statistics were meta-analysed using a standard error-weighted approach. Subsequently, we functionally annotated GWAS variants that reached genome-wide significance (P < 5 × 10-8). We also performed a gene-based GWAS, pathway analysis and linkage disequilibrium score regression and Mendelian randomization (MR) analyses. MAIN RESULTS AND THE ROLE OF CHANCE: We identified four loci associated with AMH levels at P < 5 × 10-8: the previously reported MCM8 locus and three novel signals in or near AMH, TEX41 and CDCA7. The strongest signal was a missense variant in the AMH gene (rs10417628). Most prioritized genes at the other three identified loci were involved in cell cycle regulation. Genetic correlation analyses indicated a strong positive correlation among single nucleotide polymorphisms for AMH levels and for age at menopause (rg = 0.82, FDR = 0.003). Exploratory two-sample MR analyses did not support causal effects of AMH on breast cancer or polycystic ovary syndrome risk, but should be interpreted with caution as they may be underpowered and the validity of genetic instruments could not be extensively explored. LARGE SCALE DATA: The full AMH GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Whilst this study doubled the sample size of the most recent GWAS, the statistical power is still relatively low. As a result, we may still lack power to identify more genetic variants for AMH and to determine causal effects of AMH on, for example, breast cancer. Also, follow-up studies are needed to investigate whether the signal for the AMH gene is caused by reduced AMH detection by certain assays instead of actual lower circulating AMH levels. WIDER IMPLICATIONS OF THE FINDINGS: Genes mapped to the MCM8, TEX41 and CDCA7 loci are involved in the cell cycle and processes such as DNA replication and apoptosis. The mechanism underlying their associations with AMH may affect the size of the ovarian follicle pool. Altogether, our results provide more insight into the biology of AMH and, accordingly, the biological processes involved in ovarian ageing. STUDY FUNDING/COMPETING INTEREST(S): Nurses' Health Study and Nurses' Health Study II were supported by research grants from the National Institutes of Health (CA172726, CA186107, CA50385, CA87969, CA49449, CA67262, CA178949). The UK Medical Research Council and Wellcome (217065/Z/19/Z) and the University of Bristol provide core support for ALSPAC. This publication is the work of the listed authors, who will serve as guarantors for the contents of this article. A comprehensive list of grants funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). Funding for the collection of genotype and phenotype data used here was provided by the British Heart Foundation (SP/07/008/24066), Wellcome (WT092830M and WT08806) and UK Medical Research Council (G1001357). M.C.B., A.L.G.S. and D.A.L. work in a unit that is funded by the University of Bristol and UK Medical Research Council (MC_UU_00011/6). M.C.B.'s contribution to this work was funded by a UK Medical Research Council Skills Development Fellowship (MR/P014054/1) and D.A.L. is a National Institute of Health Research Senior Investigator (NF-0616-10102). A.L.G.S. was supported by the study of Dynamic longitudinal exposome trajectories in cardiovascular and metabolic non-communicable diseases (H2020-SC1-2019-Single-Stage-RTD, project ID 874739). The Doetinchem Cohort Study was financially supported by the Ministry of Health, Welfare and Sports of the Netherlands. The funder had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript. Ansh Labs performed the AMH measurements for the Doetinchem Cohort Study free of charge. Ansh Labs was not involved in the data analysis, interpretation or reporting, nor was it financially involved in any aspect of the study. R.M.G.V. was funded by the Honours Track of MSc Epidemiology, University Medical Center Utrecht with a grant from the Netherlands Organization for Scientific Research (NWO) (022.005.021). The Study of Women's Health Across the Nation (SWAN) has grant support from the National Institutes of Health (NIH), DHHS, through the National Institute on Aging (NIA), the National Institute of Nursing Research (NINR) and the NIH Office of Research on Women's Health (ORWH) (U01NR004061; U01AG012505, U01AG012535, U01AG012531, U01AG012539, U01AG012546, U01AG012553, U01AG012554, U01AG012495). The SWAN Genomic Analyses and SWAN Legacy have grant support from the NIA (U01AG017719). The Generations Study was funded by Breast Cancer Now and the Institute of Cancer Research (ICR). The ICR acknowledges NHS funding to the NIHR Biomedical Research Centre. The content of this manuscript is solely the responsibility of the authors and does not necessarily represent official views of the funders. The Sister Study was funded by the Intramural Research Program of the National Institutes of Health (NIH), National Institute of Environmental Health Sciences (Z01-ES044005 to D.P.S.); the AMH assays were supported by the Avon Foundation (02-2012-065 to H.B. Nichols and D.P.S.). The breast cancer genome-wide association analyses were supported by the Government of Canada through Genome Canada and the Canadian Institutes of Health Research, the 'Ministère de l'Économie, de la Science et de l'Innovation du Québec' through Genome Québec and grant PSR-SIIRI-701, The National Institutes of Health (U19 CA148065, X01HG007492), Cancer Research UK (C1287/A10118, C1287/A16563, C1287/A10710) and The European Union (HEALTH-F2-2009-223175 and H2020 633784 and 634935). All studies and funders are listed in Michailidou et al. (Nature, 2017). F.J.M.B. has received fees and grant support from Merck Serono and Ferring BV. D.A.L. has received financial support from several national and international government and charitable funders as well as from Medtronic Ltd and Roche Diagnostics for research that is unrelated to this study. N.S. is scientific consultant for Ansh Laboratories. The other authors declare no competing interests.

Do sex hormones confound or mediate the effect of chronotype on breast and prostate cancer? A Mendelian randomization study.

Morning-preference chronotype has been found to be protective against breast and prostate cancer. Sex hormones have been implicated in relation to chronotype and the development of both cancers. This study aimed to assess whether sex hormones confound or mediate the effect of chronotype on breast and prostate cancer using a Mendelian Randomization (MR) framework. Genetic variants associated with chronotype and sex hormones (total testosterone, bioavailable testosterone, sex hormone binding globulin, and oestradiol) (p<5×10-8) were obtained from published genome-wide association studies (n≤244,207 females and n≤205,527 males). These variants were used to investigate causal relationships with breast (nCases/nControls = 133,384/113,789) and prostate (nCases/nControls = 79,148/61,106) cancer using univariable, bidirectional and multivariable MR. In females, we found evidence for: I) Reduced risk of breast cancer per category increase in morning-preference (OR = 0.93, 95% CI:0. 88, 1.00); II) Increased risk of breast cancer per SD increase in bioavailable testosterone (OR = 1.10, 95% CI: 1.01, 1.19) and total testosterone (OR = 1.15, 95% CI:1.07, 1.23); III) Bidirectional effects between morning-preference and both bioavailable and total testosterone (e.g. mean SD difference in bioavailable testosterone = -0.08, 95% CI:-0.12, -0.05 per category increase in morning-preference vs difference in morning-preference category = -0.04, 95% CI: -0.08, 0.00 per SD increase in bioavailable testosterone). In males, we found evidence for: I) Reduced risk of prostate cancer per category increase in morning-preference (OR = 0.90, 95% CI: 0.83, 0.97) and II) Increased risk of prostate cancer per SD increase in bioavailable testosterone (OR = 1.22, 95% CI: 1.08, 1.37). No bidirectional effects were found between morning-preference and testosterone in males. While testosterone levels were causally implicated with both chronotype and cancer, there was inconsistent evidence for testosterone as a mediator of the relationship. The protective effect of morning-preference on both breast and prostate cancer is clinically interesting, although it may be difficult to effectively modify chronotype. Further studies are needed to investigate other potentially modifiable intermediates.

Genetic insights into biological mechanisms governing human ovarian ageing.

Reproductive longevity is essential for fertility and influences healthy ageing in women1,2, but insights into its underlying biological mechanisms and treatments to preserve it are limited. Here we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in about 200,000 women of European ancestry. These common alleles were associated with clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations3. The identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR-associated genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increases fertility and extends reproductive life in mice. Causal inference analyses using the identified genetic variants indicate that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases the risk of hormone-sensitive cancers. These findings provide insight into the mechanisms that govern ovarian ageing, when they act, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

Genome-wide association study meta-analysis identifies three novel loci for circulating anti-Müllerian hormone levels in women.

Anti-Müllerian hormone (AMH) is expressed by antral stage ovarian follicles in women. Consequently, circulating AMH levels are detectable until menopause. Variation in age-specific AMH levels has been associated with breast cancer and polycystic ovary syndrome (PCOS), amongst other diseases. Identification of genetic variants underlying variation in AMH levels could provide clues about the physiological mechanisms that explain these AMH-disease associations. To date, only one variant in MCM8 has been identified to be associated with circulating AMH levels in women. We aimed to identify additional variants for AMH through a GWAS meta-analysis including data from 7049 premenopausal women of European ancestry, which more than doubles the sample size of the largest previous GWAS. We identified four loci associated with AMH levels at p < 5×10 -8 : the previously reported MCM8 locus and three novel signals in or near AMH, TEX41 , and CDCA7 . The strongest signal was a missense variant in the AMH gene (rs10417628). Most prioritized genes at the other three identified loci were involved in cell cycle regulation. Genetic correlation analyses indicated a strong positive correlation among SNPs for AMH levels and for age at menopause (r g = 0.82, FDR=0.003). Exploratory Mendelian randomization analyses did not support a causal effect of AMH on breast cancer or PCOS risk, but should be interpreted with caution as they may be underpowered and the validity of genetic instruments could not be extensively explored. In conclusion, we identified a variant in the AMH gene and three other loci that may affect circulating AMH levels in women.

Large Copy-Number Variants in UK Biobank Caused by Clonal Hematopoiesis May Confound Penetrance Estimates.

Large copy-number variants (CNVs) are strongly associated with both developmental delay and cancer, but the type of disease depends strongly on when and where the mutation occurred, i.e., germline versus somatic. We used microarray data from UK Biobank to investigate the prevalence and penetrance of large autosomal CNVs and chromosomal aneuploidies using a standard CNV detection algorithm not designed for detecting mosaic variants. We found 160 individuals that carry >10 Mb copy number changes, including 56 with whole chromosome aneuploidies. Nineteen (12%) individuals had a diagnosis of Down syndrome or other developmental disorder, while 84 (52.5%) individuals had a diagnosis of hematological malignancies or chronic myeloproliferative disorders. Notably, there was no evidence of mosaicism in the blood for many of these large CNVs, so they could easily be mistaken for germline alleles even when caused by somatic mutations. We therefore suggest that somatic mutations associated with blood cancers may result in false estimates of rare variant penetrance from population biobanks.

Using human genetics to understand the disease impacts of testosterone in men and women.

Testosterone supplementation is commonly used for its effects on sexual function, bone health and body composition, yet its effects on disease outcomes are unknown. To better understand this, we identified genetic determinants of testosterone levels and related sex hormone traits in 425,097 UK Biobank study participants. Using 2,571 genome-wide significant associations, we demonstrate that the genetic determinants of testosterone levels are substantially different between sexes and that genetically higher testosterone is harmful for metabolic diseases in women but beneficial in men. For example, a genetically determined 1 s.d. higher testosterone increases the risks of type 2 diabetes (odds ratio (OR) = 1.37 (95% confidence interval (95% CI): 1.22-1.53)) and polycystic ovary syndrome (OR = 1.51 (95% CI: 1.33-1.72)) in women, but reduces type 2 diabetes risk in men (OR = 0.86 (95% CI: 0.76-0.98)). We also show adverse effects of higher testosterone on breast and endometrial cancers in women and prostate cancer in men. Our findings provide insights into the disease impacts of testosterone and highlight the importance of sex-specific genetic analyses.

Genetic studies of accelerometer-based sleep measures yield new insights into human sleep behaviour.

Sleep is an essential human function but its regulation is poorly understood. Using accelerometer data from 85,670 UK Biobank participants, we perform a genome-wide association study of 8 derived sleep traits representing sleep quality, quantity and timing, and validate our findings in 5,819 individuals. We identify 47 genetic associations at P < 5 × 10-8, of which 20 reach a stricter threshold of P < 8 × 10-10. These include 26 novel associations with measures of sleep quality and 10 with nocturnal sleep duration. The majority of identified variants associate with a single sleep trait, except for variants previously associated with restless legs syndrome. For sleep duration we identify a missense variant (p.Tyr727Cys) in PDE11A as the likely causal variant. As a group, sleep quality loci are enriched for serotonin processing genes. Although accelerometer-derived measures of sleep are imperfect and may be affected by restless legs syndrome, these findings provide new biological insights into sleep compared to previous efforts based on self-report sleep measures.

Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms.

Being a morning person is a behavioural indicator of a person's underlying circadian rhythm. Using genome-wide data from 697,828 UK Biobank and 23andMe participants we increase the number of genetic loci associated with being a morning person from 24 to 351. Using data from 85,760 individuals with activity-monitor derived measures of sleep timing we find that the chronotype loci associate with sleep timing: the mean sleep timing of the 5% of individuals carrying the most morningness alleles is 25 min earlier than the 5% carrying the fewest. The loci are enriched for genes involved in circadian regulation, cAMP, glutamate and insulin signalling pathways, and those expressed in the retina, hindbrain, hypothalamus, and pituitary. Using Mendelian Randomisation, we show that being a morning person is causally associated with better mental health but does not affect BMI or risk of Type 2 diabetes. This study offers insights into circadian biology and its links to disease in humans.

Genome-wide association study of anti-Müllerian hormone levels in pre-menopausal women of late reproductive age and relationship with genetic determinants of reproductive lifespan.

Anti-Müllerian hormone (AMH) is required for sexual differentiation in the fetus, and in adult females AMH is produced by growing ovarian follicles. Consequently, AMH levels are correlated with ovarian reserve, declining towards menopause when the oocyte pool is exhausted. A previous genome-wide association study identified three genetic variants in and around the AMH gene that explained 25% of variation in AMH levels in adolescent males but did not identify any genetic associations reaching genome-wide significance in adolescent females. To explore the role of genetic variation in determining AMH levels in women of late reproductive age, we carried out a genome-wide meta-analysis in 3344 pre-menopausal women from five cohorts (median age 44-48 years at blood draw). A single genetic variant, rs16991615, previously associated with age at menopause, reached genome-wide significance at P = 3.48 × 10-10, with a per allele difference in age-adjusted inverse normal AMH of 0.26 standard deviations (SD) (95% confidence interval (CI) [0.18,0.34]). We investigated whether genetic determinants of female reproductive lifespan were more generally associated with pre-menopausal AMH levels. Genetically-predicted age at menarche had no robust association but genetically-predicted age at menopause was associated with lower AMH levels by 0.18 SD (95% CI [0.14,0.21]) in age-adjusted inverse normal AMH per one-year earlier age at menopause. Our findings provide genetic support for the well-established use of AMH as a marker of ovarian reserve.

Mosaic Turner syndrome shows reduced penetrance in an adult population study.

PURPOSE: Many women with X chromosome aneuploidy undergo lifetime clinical monitoring for possible complications. However, ascertainment of cases in the clinic may mean that the penetrance has been overestimated. METHODS: We characterized the prevalence and phenotypic consequences of X chromosome aneuploidy in a population of 244,848 women over 40 years of age from UK Biobank, using single-nucleotide polymorphism (SNP) array data. RESULTS: We detected 30 women with 45,X; 186 with mosaic 45,X/46,XX; and 110 with 47,XXX. The prevalence of nonmosaic 45,X (12/100,000) and 47,XXX (45/100,000) was lower than expected, but was higher for mosaic 45,X/46,XX (76/100,000). The characteristics of women with 45,X were consistent with the characteristics of a clinically recognized Turner syndrome phenotype, including short stature and primary amenorrhea. In contrast, women with mosaic 45,X/46,XX were less short, had a normal reproductive lifespan and birth rate, and no reported cardiovascular complications. The phenotype of women with 47,XXX included taller stature (5.3 cm; SD = 5.52 cm; P = 5.8 × 10-20) and earlier menopause age (5.12 years; SD = 5.1 years; P = 1.2 × 10-14). CONCLUSION: Our results suggest that the clinical management of women with 45,X/46,XX mosaicism should be minimal, particularly those identified incidentally.

Using genetics to understand the causal influence of higher BMI on depression.

BACKGROUND: Depression is more common in obese than non-obese individuals, especially in women, but the causal relationship between obesity and depression is complex and uncertain. Previous studies have used genetic variants associated with BMI to provide evidence that higher body mass index (BMI) causes depression, but have not tested whether this relationship is driven by the metabolic consequences of BMI nor for differences between men and women. METHODS: We performed a Mendelian randomization study using 48 791 individuals with depression and 291 995 controls in the UK Biobank, to test for causal effects of higher BMI on depression (defined using self-report and Hospital Episode data). We used two genetic instruments, both representing higher BMI, but one with and one without its adverse metabolic consequences, in an attempt to 'uncouple' the psychological component of obesity from the metabolic consequences. We further tested causal relationships in men and women separately, and using subsets of BMI variants from known physiological pathways. RESULTS: Higher BMI was strongly associated with higher odds of depression, especially in women. Mendelian randomization provided evidence that higher BMI partly causes depression. Using a 73-variant BMI genetic risk score, a genetically determined one standard deviation (1 SD) higher BMI (4.9 kg/m2) was associated with higher odds of depression in all individuals [odds ratio (OR): 1.18, 95% confidence interval (CI): 1.09, 1.28, P = 0.00007) and women only (OR: 1.24, 95% CI: 1.11, 1.39, P = 0.0001). Meta-analysis with 45 591 depression cases and 97 647 controls from the Psychiatric Genomics Consortium (PGC) strengthened the statistical confidence of the findings in all individuals. Similar effect size estimates were obtained using different Mendelian randomization methods, although not all reached P < 0.05. Using a metabolically favourable adiposity genetic risk score, and meta-analysing data from the UK biobank and PGC, a genetically determined 1 SD higher BMI (4.9 kg/m2) was associated with higher odds of depression in all individuals (OR: 1.26, 95% CI: 1.06, 1.50], P = 0.010), but with weaker statistical confidence. CONCLUSIONS: Higher BMI, with and without its adverse metabolic consequences, is likely to have a causal role in determining the likelihood of an individual developing depression.

Genome-Wide and Abdominal MRI Data Provide Evidence That a Genetically Determined Favorable Adiposity Phenotype Is Characterized by Lower Ectopic Liver Fat and Lower Risk of Type 2 Diabetes, Heart Disease, and Hypertension.

Recent genetic studies have identified alleles associated with opposite effects on adiposity and risk of type 2 diabetes. We aimed to identify more of these variants and test the hypothesis that such favorable adiposity alleles are associated with higher subcutaneous fat and lower ectopic fat. We combined MRI data with genome-wide association studies of body fat percentage (%) and metabolic traits. We report 14 alleles, including 7 newly characterized alleles, associated with higher adiposity but a favorable metabolic profile. Consistent with previous studies, individuals carrying more favorable adiposity alleles had higher body fat % and higher BMI but lower risk of type 2 diabetes, heart disease, and hypertension. These individuals also had higher subcutaneous fat but lower liver fat and a lower visceral-to-subcutaneous adipose tissue ratio. Individual alleles associated with higher body fat % but lower liver fat and lower risk of type 2 diabetes included those in PPARG, GRB14, and IRS1, whereas the allele in ANKRD55 was paradoxically associated with higher visceral fat but lower risk of type 2 diabetes. Most identified favorable adiposity alleles are associated with higher subcutaneous and lower liver fat, a mechanism consistent with the beneficial effects of storing excess triglycerides in metabolically low-risk depots.